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I.  INTRODUCTION 

The  horizontal  plane  stability  of  towed  vessels  is  a  very  important  field.  If 
the  motion  of  a  vessel  under  tow  is  unstable,  a  very  dangerous  situation  exists; 
the  motion  of  a  several  thousand  ton  vessel  could  endanger  the  ship  towing  it 
[Refs.  1,  2,  3].  A  method  needs  to  be  developed  not  only  to  ensure  the  stability  of 
the  towed  vessel's  motion,  but  to  control  it  [Refs.  4,  5].  If  the  motion  of  the  towed 
vessel  can  be  controlled,  the  next  step  is  to  optimize  its  control,  both  for  safety 
and  fuel  considerations.  Usually  the  only  means  of  control  of  a  vessel  under  tow 
is  through  the  towline.  Under  these  conditions,  stability  of  the  vessel  is  marginal; 
in  some  cases,  the  vessel  is  actually  unstable. 

This  thesis  will  develop  the  equations  of  motion  for  a  vessel  in  tow.  It  will 
also  state  and  describe  the  unique  forces  associated  with  a  vessel  being  towed. 
Finally,  it  will  describe  how  state-space  theory  is  being  applied  to  control  the 
towed  vessel's  motion  through  the  movement  of  the  attachment  point  of  the 
towline  along  the  beam. 

A  control  system  based  on  sliding  mode  control  and  the  linear  quadratic 
regulator  for  the  linearized  equations  of  motion  will  be  presented,  with  results  in 
non-dimensional  form.  The  same  control  system  for  the  non-linearized  equations 
of  motion  will  also  be  presented,  along  with  the  results.  Comparison  of  results  for 
the  linearized  and  non-linearized  equations  of  motion  will  help  determine  the 


robustness  of  the  system.  Matrix,  will  be  used  to  calculate  the  required  gains  for 
the  system,  and  the  simulation  will  be  done  using  a  Fortran  program. 


II.  EQUATIONS  OF  MOTION 

The  two  equations  of  motion  will  be  developed  using  a  body-fixed  frame  of 
reference.  The  origin  of  this  reference  frame  is  at  the  center  of  gravity  of  the 
towed  vessel.  In  general,  the  vessel  has  six  degrees  of  freedom: 

1.  surge  along  the  xg-axis 

2.  roll  about  the  x -axis 

3.  sway  along  the  yg-axis 

4.  pitching  about  the  yg-axis 

5.  heave  along  the  zg-axis 

6.  yaw  about  the  z  -axis 

This  thesis  will  only  deal  with  motion  in  the  horizontal  plane.  Also,  a  constant 
surge  velocity  along  the  xg-axis  will  be  assumed.  Other  assumptions  include: 

1.  no  towed  vessel /towing  vessel  interaction 

2.  massless,  inextensible  towline 

3.  no  wind,  wave,  or  current  disturbances 

4.  the  body  is  symmetric 


Since  motion  is  confined  to  the  horizontal  plane,  the  applicable  equations  are  the 
sway  and  yaw  equations.  These  equations  are  standard  Principles  of  Naval 
Architecture,  or  PNA,  equations.  They  are  as  follows: 
Sway: 


m[v +  ur +wp +x  (pq+f)-y  (p2  +r2)  +z  (qr-p)] 
=  (W-B)  cos6  cos<}>  +  Y 


(la) 


Yaw: 


m[x  (v+ur+w)-y  (u-vr+wcj)]  uW 

=  ( x  w  -  x.B  )  cosG  cos<t)  +(y  W  -  ypB  )  sin0  +N 

If  we  make  the  above  assumptions  in  terms  of  the  variables  in  the  equations: 

1.  No  vertical  motion;  w=0 

2.  W=B 

3.  Small  angle  motion  (linearized  simulation  only) 

4.  No  pitch  rate;  also  negligible  roll;  0=0,  6  =  0,  <j>=0 

5.  Forward  speed  equals  nominal  speed;  u=U 

6.  Center  of  Buoyancy=Center  of  gravity;  xg,  yg=Xp,  yp 


Then  the  sway  and  yaw  equations  become: 


m(v  +  ur  +x  f)  =  Y,  (2a) 

if 


I  f  +mx  (v  +  ur)  =  N,  (2b) 

z  g  f 


The  hydrostatic  forces  on  the  right  hand  side  of  the  equations  are  then  expanded 
in  a  standard  Taylor  series: 


Y  =  Yv  +  Yr  +  Y.v  +  Y.r  +  Y    .  (3a) 

t  v  r  v  r  touxine 


N=Nv+Nr+N.v+N.f  +  N,  ..  (3b) 

F  v  r  v  r  tovAme 


where        Yv  =  J*I  |  ^c  ,   etc 
ov 


Note  that  the  rudder  terms,  N6  and  Y6  have  been  omitted,  as  the  towed  vessel's 
rudder  is  amidships;  there  is  no  applied  force  due  to  the  rudder.  The  only  applied 
force  is  due  to  the  towline.  Substituting  (3a)  into  (2a)  and  (3b)  into  (2b)  and 
rearranging  the  terms  yields: 
Sway: 

(m-Y  )v  +  (-Y.)f  =  Yv  +  (Y  -mu)r  +  Y,  ..  (4a) 

v  r  v  r  tcruAmt 

Yaw: 

(I  -N.)f-N.v=Nv+Nr+Nt   ..  (4b) 

z  r  v  v  r  towime 


which  can  then  be  rearranged  into  working  form  as 

Y  v  +  ( Y.  -  m  yb  +  ( Y  -  mu  )ty  +  YJt  =  Y   ..  (5a) 

Nv+N.ti+Ntf  +  iN. '-I  )iif  =  N   ,.  (5b) 

where  r^,  f=\\f.  All  that  remains  is  to  describe  the  towline  forces  in  greater 
detail. 

A.      TOWLINE  FORCE 

In  Figure  \,  P  is  the  attachment  point  of  the  tow.  Xp  and  yp  denote  the 
distance  from  the  towed  vessel's  center  of  gravity,  again  using  the  standard  PNA 
body-fixed  coordinate  system.  From  the  diagram,  we  can  see  that  the  horizontal 
force  associated  with  the  towline  in  the  sway  direction  is: 

T  sin  ( y  +  \|/ ) 
The  moment  causing  the  vessel  to  yaw  is: 

-Tx  ( sin  (y  +  y ) )  -  Ty  cos  ( y  +  \\f ) 
Assuming  y  and  y  are  small  angles: 

y  +x  \|/  +  y 
sin(y  +  y)=Y  +  \}/  = _ +  y 


f 

x 

1 

+ 

V 

L 

V 

J 

y    ^v 
=  \i/  i+_i  +_  +  _ 

L       L 


Path 


Also, 


Figure  1.  Towed  Vessel  Coordinate  Description. 


y  -v  +  wsiny  =  v  +u\j/ 


z;  =  y  -w\j/ ;  i>  =  y  -w\|/ 


Finally,  substituting  for  v,   in  equations  (5a)  and  (5b)  and  rearranging  terms: 
Sway: 


Yaw: 


(Yv  -m)y  +YJ,  -  -y  +  Y„y  -(Ym  -Y^)y 


/ 


-[Yu  +T 


V 


x 

1  +-1 
L 


¥  = 


Ty. 


T 


NJJ+NJ~  jV  +(N^  -Jz)Vjif-(N.M  -Ny)\j/ 


-  [  N  u  +  7> 


1  +.1 


lv 


(6a) 


(6b) 


=  —  a:  v  +  Ty 


These  two  equations  are  motion  describe  the  towed  vessel's  motion.  Note  that  the 
coordinate  y  describes  the  lateral  distance  from  the  towing  path;  yp  describes  the 
distance  of  the  towline  attachment  point  from  the  centerline  of  the  towing  vessel. 
It  is  for  yp  that  the  control  law  is  being  developed. 
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B.       ARRANGEMENT  INTO  STATE-SPACE 

Let  X]=y,  x2=\j/,  x3=y,  x4=\|/,  then,  substituting  into  equations  (6a)  and  6b) 


(Y.-m)x+Y.x. 

v  3  v     4 


(  \ 

X 

L 

J 


=  -x+[Yu+T   \+JL    }x-Yx+(Y.u-Y)x+-y 

tIp  t         j      2  v    3  v  t|r         4         r     "p 

V  J 


(7a) 


N  x*  +  (N-l  )xk 

v    3  y         z  '      4 


T 

-jcr  +    Nw+Tx 

^    p  i  p  P 


X-  -N  x~+(N.u  -N  )x+  —  x  y  +Ty 

2  p3p  ty'4         r       p^p  •?;; 


From  the  definition  of  the  state  variables,  it  can  be  seen  that 


(7b) 


x  -  x 

Al        A3 


(8a) 


and 


x  =  x 

X2       A4 


(8b) 


These  comprise  the  first  two  of  four  state  equations.  The  second  two  state 
equations  must  come  from  equations  (7a)  and  (7b)  above.  Note  that  these 
equations  are  coupled  in  x3  and  x4.  Algebraic  elimination  yields  the  final  two 
equations  of  motion.  Multiplying  equation  (7a)  by  (N  -lz)  and  equation  (7b) 
by   (  -  Yy  )   yields  an  equation  for  x3 . 


[(N.-l  )(Y.-m)-N.y.]x. 
=  Z[(N.  -/  )-xYJx, 


y  m  +t 


f        \ 

X 

1  +  _£ 

L 


\ 


JJ 


(N.  -J  ) 

V  2 


-y. 


N  u  +Tx 


(  \ 

x 
1  +  _! 


+  [  n  y .  -  y  (  n.  - 1 )  ]  *, 

+  [  (  y.m  -  y,  (  n.  -  /  )  -  y.  (  n.m  -  n,  ]  *4 


(N.  -/  ) 

¥     z    -  y 

r  * 


(  \ 

x 

1  +_£ 


(8c) 


Similarly,  multiplying  equation  (7a)  by    (  -Nu )   and  equation  (7b)  by  (  Y9  -  m ) 
and  adding  yields  an  equation  for  x4 . 
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[(N.  -/  )(Y.  -m)-N.yJxA 


T 


_[*  (Y.  -m)-N.]x, 

{(Y  -m) 


N  u  +Tx 

v  p 


(         \ 

X 

L 

v  yj 


-N. 


Yw  +T 


a: 

L 

v  yj 


}*, 


+  [N  Y  -N  (Y.  -m)]x, 


V         V  V  V 


+  [  -N  mu  +N,  m  -N,  Y.  +N  y,  )*, 


+  T 


(Y-m) 


(         > 

X 

N. 

1  +  jl 

V 

L 

~~L 

v           y 

-r  y* 


(8d) 


Equations  (8a),  (8b),  (8c)  and  (8d)  can  also  be  expressed  in  state-space  matrix  form 


p- 

*1 

*2 

= 

*3 

.*«_ 

0      0      10 
0      0      0      1 


fl       fl       fl       fl 

31         32        33        34 


fl41      fl42     fl43     *44 


r      -. 

~ 

*1 

0 

X2 

0 

X3 

+ 

K 

lX*. 

K_ 

Where, 


D  =  (  Ny  -  Iz  )    (  Y,  -  m  )   -  N,  Yv 


a3L=D-1-I  [    (Nv-Iz)   -xpYv 


(9) 


A32=D"M 


y,  ,u  +  r 


l  +  _p 
v         L   J 


(  N„  -  I J   -Yw 


Nu  +  Txr 


1  +  _£ 
V  L    J 
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a3,=D-:  [NvYy-Yv  (%-  Iz)   ] 

a5<  =  D--  [    (Yvu-Yy)    (Ny-Iz)   -Y^(N^u-N^)   ] 


a,.=D---±[x(Y-m)-N„] 


L 


C      v  V 


a42  =£)--{  (  Yv-m) 


N..u  +  Tx. 


1  +_£ 

K  L    , 


-  N... 


y„u  +  r 


(     Ol 

1  +  -E. 

I    L  )\ 


a43  =  ZT:  [  NVYV-NV  (  Yv-m)   ] 


a«  =  D--  [  NvYy  -  NyYv  +  (N^-Nvu)  m] 


b,  =  d~:  r 


(V^J.y 


1    +    _£ 

v         L  > 


£4  =  D-'-T 


(  y*  -  m  ) 


(    o 

1  +_z 

V    L  J 


Hi 

L 
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III.  CONTROL  SYSTEM  DESIGN 

A.      SYSTEM  AUGMENTATION 

Equation  (9)  defines  the  state-space  form  of  the  equations  of  motion.  One 
more  state  must  be  defined  and  added  to  this  equation;  that  of  the  ordered,  or 
commanded,  yp.  This  variable  shall  be  called  y^..  The  usual  response  of  yp  with 
respect  to  time  should  look  like  a  first  order  system,  as  in  Figure  2: 


'pc 


0 


-.5 


time 


Figure  2.  y   Response  with  Respect  to  Time. 
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The  equation  governing  this  response  is  as  follows: 

-1        1 

or     v  =  —  y  +  —  v 
p  p 

where  Tp  is  the  non-dimensional  time  constant  for  the  towline  attachment  point 
motion  control  system.  If  y^  is  now  considered  the  input  to  the  system,  and  y_ 
is  considered  a  state  variable,  the  matrix  equation  becomes: 


-1 

Y 


0      0     0     0 


r       _ 

y? 

"l    " 

Y 

x, 

V 

1 

0 

X2 

+ 

0 

X, 

3 

0 

l.V 

0 

(10) 


P<= 


0  0      0      10 

0  0      0      0      1 

h  a  a  a      a 

u3  31        32        33        34 

b  a  a  a     a 

4  41         42        43        44 


where  a31,  a32,  a^,  a^,  a41,  a^,  a43,  a4,  b3  and  b4  are  as  previously  defined.  In  a 
sense,  the  actual  distance  of  the  towline  attachment  point  yp  is  treated  as  an  extra 
state  to  the  system.  The  control  system  will  be  designed  based  on  this  augmented 
system  of  equations.  The  control  system  will  be  designed  based  on  two  principles 
of  state-space  theory;  Linear  Quadratic  Regulator,  and  Sliding  Mode  Control. 

B.       LINEAR  QUADRATIC  REGULATOR 

The  linear  quadratic  regulator  function  arose  from  efforts  to  find  an 
optimum  means  of  control  [Ref.  6].  Normally,  a  gain  matrix  for  a  particular 
control  system  is  calculated  to  achieve  specific  closed-loop  pole  locations.  With 
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the  linear  quadratic  regulator,  or  LQR,  a  specific  performance  criterion  J  (or  "cost 
function")  is  defined,  with  the  only  stipulations  on  the  poles  being  that  they  be 
negative,  or  stable.  This  criterion  J  is  defined  as: 

J  =  j[x'(t)Q(t)x(t)+u,(t)R(t)u(t)]dt 
The  matrices  Q  and  R  are  weighting  matrices.  The  matrix  Q  is  a  state 
weighting  matrix,  and  R  is  a  control  weighting  matrix.  The  gains  for  a  control 
system  can  be  calculated  based  on  the  defined  matrices  Q  and  R.  For  instance,  if 
the  elements  of  the  matrix  Q  are  small  relative  to  R,  the  system  will  tolerate  large 
errors  in  the  final  state  with  very  little  control  effort.  Conversely,  if  Q  is  large 
relative  to  R,  very  small  errors  in  the  state  X  will  result,  but  with  considerable 
control  effort. 

C      SLIDING  MODE  CONTROL 

The  second  aspect  of  state-space  theory  to  be  used  in  controlling  the  towed 
vessel  is  that  of  sliding  mode  control.  Since  the  control  law  is  based  on  a 
linearized  set  of  equations  of  motion,  a  lot  of  uncertainty  in  the  response  exists. 
Also,  some  of  the  parameters  may  vary.  A  good  example  for  the  towed  vessel 
system  is  the  towline  tension  T,  which  will  certainly  vary  with  time.  A  control 
law  needs  to  be  developed  that  will  take  into  account  both  the  uncertainties  in  the 
parameters,  and  any  dynamics  that  have  either  not  been  modeled,  or  that  have 
been  linearized.  A  sliding  mode  control  law  can  ensure  both  stability  and 
robustness  of  the  system,  with  the  emphasis  on  robustness.  The  LQR  gains  will 
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ensure  stability  of  the  system.  Sliding  mode  control  is  ideally  suited  to  a  system 
where  the  response  oscillates  between  set  values,  such  as  the  motion  of  a  towed 
vessel;  it  uses  a  high  speed  "toggling"  control  law  to  drive  the  system  onto  a 
desired  "sliding  plane." 

Sliding  mode  control  takes  the  standard  state-space  system 

x  =  [A]x  +  [b]u  (12) 

and  defines  a  sliding  plane 

8(x)  =s1x1  +s2x2  +s3x3  +s4x4  (13) 

and  the  coefficient  Sj  is  arbitrary.  Equation  (13  )  can  be  written  as 

sTx  =  0 
where 

S       —  LSj/^2'    3'    4  J 

Next,  define  the  Lyapunov  function 

V(x)=i[e(x)]2 
2 

stability  is  guaranteed,  provided  v  (x)  is  a  negative  definite  function.  Another 
way  to  express  this  is 

V(x)=Q(x)S(x)  =  -t\2Q(x)  (14) 
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therefore, 

6  =  -r|2  sign  6  (15) 

Since 

Q(x)=sT(x),b(x)=sT(x)=sT(Ax+bu) 
so,  substituting  for  0  in  Equation  (15): 

s  T  ( Ax  +  bu )  =  -r|2  sign  (0 ) 
and  solving  for  u: 

u  =  -(sTb)-' sTAx  -r\2 (sTbY' sign(Q)  <16> 

Equation  (16)  is  a  sliding  mode  control  law  for  a  generic  system. 

D.      THE  CONTROL  LAW 

All  that  remains  is  to  define  a  control  law  for  the  augmented  equations  of 
motion  for  the  towed  vessel.  The  control  law  is  denned  using  sliding  mode 
control  and  linear  quadratic  regulator  principles.  Take  the  augmented  system: 
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(17) 


and  split  the  augmented  matrix  A  into  Au,  A12,  A21/  and  A^  as  shown.  Also  split 
the  state  vector  into 


yi=y„;y2 -[xlfx2,Xyx4] 


Rewrite  equation  (17)  as 


y2  =  A2^+A2iy2  +  b2u 

The  original  augmented  state  vector  x  now  equals  (y^  y2). 
Next,  define  the  sliding  plane  for  the  towed  vessel  as 


eU)=e(yi/i/2)=i/1+s'y2=0 


and 


6  =  -n2s/g«(e)  =  i/1+sTi/2 


(18) 


(19) 
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but 


yi=Anyi+fciw  and  y2=A2iyi+A22y2 

so,  substituting  for  y.  and  y2  is  equation  (19): 


Any1  +  &,«  +sT(A21ya  *A22y2)  =  -n2si'g«(e) 


solving  for  u: 


1 
-l(Au  +sTAn)y^  ^sTA12y2]  -      j]2sign(Q) 


-1 

T 


i 


(20) 


This  is  the  system  control  law  for  the  towed  vessel. 

The  gains  [s]  for  the  sliding  mode  control  law  from  minimization  of  the 
linear  quadratic  regulator  cost  function  J.  Recall  that 


J  =  j{[xY[Q][x}+[u]'[R][u])dt 
in  Equation  (21)  the  state  weighting  matrix  Q  is  as  follows: 


(21) 


Q  = 


0 

0 

0    0" 

0 

<?22 

0    0 

0 

0 

0    0 

0 

0 

0   0_ 

where  q^  represents  an  error  in  radians  from  the  path  of  the  towing  vessel.  We 
are  choosing  to  weigh  the  state  variable  x2  =  \y.  The  control  weighting  matrix  R 
will  be  defined  by 
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K  =  (ot)-2 

where  a,  the  maximum  non-dimensional  distance  yp  will  be  weighted.  With  Q 
and  R  defined,  the  cost  matrix  function  becomes: 


J^minjiq^+Ru^dt  (22) 


or 


/  =  min  f  ( qny*  +  Ry2p  ) dt  <23) 

The  smaller  the  value  for  R  is,  the  smaller  the  control  effort  required,  but  a  larger 
state  error  in  q^  will  be  have  to  be  accepted  as  a  trade-off.  Conversely,  the  larger 
the  value  for  R  is  the  more  control  required,  with  very  little  state  error  in  q^. 

The  resulting  gains  are  then  placed  into  the  control  law  for  the  augmented 
system  of  equations,  Equation  (20),  and  the  response  of  the  system  is  obtained. 
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IV.  RESULTS 

A.  MATRIXX  AND  FORTRAN  IMPLEMENTATION 

MatrixX  is  an  outstanding  tool  for  understanding  the  response  of  a  control 
system.  A  Fortran  program  was  written  to  use  with  MatrixX.  Since  the  coefficients 
of  the  [A]  and  [B]  matrices  of  the  state-space  equation  are  constant  with  time, 
they  can  be  put  directly  into  the  MatrixX  software,  and  the  resulting  Linear 
Quadratic  Regulator  gains  are  easily  calculated.  Without  the  use  of  MatrixX, 
solving  the  linear  quadratic  integral  is  a  very  difficult  problem  in  numerical 
analysis.  The  actual  simulation  of  the  system  response  was  accomplished  with  a 
Fortran  program  on  the  VAX  computer.  The  LQR  gains  obtained  from  MatrixX 
are  put  into  a  data  file,  along  with  the  non-dimensional  parameters  of  Tp,  LL,  Xp, 
T,  and  the  initial  conditions  of  Y,  \|/,  and  the  maximum  distance  Yp  can  travel 
athwartships.  The  simulation  then  reads  this  data  file,  computes  the  sliding  mode 
control  and  plots  the  response  of  the  parameters  Y,  Yp,  y,  and  Y^  with  respect  to 
time. 

B.  INITIAL  CONDITIONS  AND  CONSTANT  PARAMETERS 

The  results  are  presented  graphically,  and  are  in  non-dimensional  form,  with 
the  exception  of  the  yaw  angle  \y,  which  is  in  degrees.  The  non-dimensionalization 
is  standard  Principles  of  Naval  Architecture:  the  relationships  describing  the 
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nondimensional  parameters  are  included  in  the  Appendix.  Two  vessels  were 
studied: 

1.  528  foot  mariner 

2.  1066.3  foot  tanker 

A  summary  of  the  hull  particulars  and  hydrodynamic  derivatives  for  both  these 
vessels  in  also  included  in  the  Appendix.  Four  knots  was  used  as  the  nominal 
forward  velocity.  The  mariner  is  stable  at  four  knots;  the  tanker  is  unstable.  The 
towline  tension  for  both  vessels  is  taken  from  resistance  curves  at  four  knots,  and 
then  non-dimensionalized;  the  value  used  is  0.001  Xp,  the  longitudinal  distance 
from  the  towed  vessel  center  of  gravity  to  the  towline  attachment  point,  is 
assumed  to  be  constant.  A  value  of  0.5  is  used.  The  time  required  for  Yp  to 
"match"  Y^  as  in  Figure  2,  T  is  set  as  0.5.  The  maximum  non-dimensional 
distance  Yp  can  travel  from  port  to  starboard  extremes  is  taken  to  be  0.1;  therefore 
the  maximum  distance  from  an  extreme  to  centerline  is  0.05.  The  initial  conditions 
for  each  data  run  are  Y=0.5,  y=0,  and  Yp=0.05  as  in  Figure  3.  Also,  the  weighting 
matrices  were  set  at  q22=5°  (or  .087  radians)  for  Q,  the  state  error  weighting 
matrix,  and  a=0.015  for  R,  the  control  effect  weighting  for  Y^..  Four  plots  were 
generated  for  each  simulation.  The  y-distance  from  the  path  (Y),  the  towline 
attachment  point  lateral  offset  (Yp),  the  towline  attachment  point  offset  command 
(Y^)  and  the  yaw  angle  (y)  are  all  plotted  vs  nondimensional  time.  Variance  of 
parameters  of  interest  are  shown  on  subsequent  runs. 


22 


Towed 
Vessel 


Towing 
Vessel 


Figure  3.  Initial  Conditions. 
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1.  Figures  4  and  5:  Standard  Mariner  and  Tanker 

Figure  4  shows  the  response  of  the  mariner  with  the  "standard" 
parameters  and  initial  conditions  stated  above;  Figure  5  shows  the  "standard" 
tanker.  Overshoot  in  Y  and  y  exists  for  both  cases,  although  it  is  higher  in  the 
mariner.  The  Y^  graph  shows  that  significantly  more  effort  is  required  to  control 
the  tanker,  and  the  oscillation  in  Yp  is  much  greater  than  for  the  mariner.  This 
makes  sense,  as  the  tanker  is  inherently  more  unstable.  Both  vessels  converge  to 
a  straightline  path  in  approximately  60  non-dimensional  seconds.  As  the  tanker 
is  more  unstable,  it  will  be  used  to  demonstrate  all  other  variations  in  parameters. 
Henceforth  all  references  to  "normal"  or  "standard"  tanker  response  will  be  to  the 
response  shown  in  Figure  5. 

2.  Figure  6:  Tanker  with  Y,  Y  Not  Observable 

Figure  6  shows  the  vessel  response  if  the  state  variables  Y  and  Y  are  not 
observable  and  are  not  used  in  the  control  law.  Both  Y  and  Y  converge  but  the 
time  required  is  significant;  over  100  seconds.  In  this  case  there  is  no  overshoot 
in  Y,  but  considerable  overshoot  in  \\f.  \\f,  Y^  and  Yp  all  oscillate  considerably;  but 
in  the  end,  the  system  does  converge. 

3.  Figures  7,  8:  Observable  Tanker,  Variance  of  6(x) 

Recall  from  Equation  (13)  that  the  sliding  plane  6(x)  is  in  part  arbitrarily 
defined,  as  the  gain  rj2  is  chosen  by  the  designer.  For  the  standard  runs,  rj2  was 
set  to  be  1.0.  Figure  7  shows  the  tanker  response  if  rj2  in  6(x)  is  set  to  0.5,  with  all 
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Figure  4.  Mariner  Response  with  Standard  Parameters. 
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Figure  5.  Tanker  Response  with  Standard  Parameters. 
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Figure  6.  Tanker  with  Y,  Y  Not  Observable. 
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state  variables  observable.  The  system  still  converges  in  adequate  time;  the 
parameters  Y,  Ypc/  and  Yp  behave  as  in  the  standard  tanker.  There  is 
approximately  twice  as  much  overshoot  in  y,  however.  In  Figure  8,  rj2  is  set  to 
2.0;  this  brought  the  overshoot  in  y  down  to  the  level  of  the  standard  observable 
tanker,  without  altering  the  favorable  response  of  the  other  three  parameters. 

4.  Figures  9,  10,  11:  Non-Observable  Tanker,  Variance  of  6(x) 

The  state  variables  Y  and  Y  were  assumed  to  be  not  observable  in 
Figures  9  and  10.  In  Figure  9,  rj2  was  set  to  0.5,  with  disastrous  results.  The  towed 
vessel  becomes  unstable.  All  four  parameters  diverge  with  time.  In  Figure  10, 
however,  rj2  in  0(x)  is  set  to  2.0,  and  the  system  response  stabilizes.  The  overshoot 
in  \\f  is  approximately  as  the  same  as  the  standard  tanker,  but  the  overshoot  in  Yp 
and  Y^  is  twice  that  of  the  normal  tanker.  The  settling  time  in  Y  is  double  that 
of  the  normal  tanker,  but  the  settling  time  in,  Y^.  and  Yp  is  cut  to  forty  seconds. 
Obviously  there  is  a  value  between  rj2  =  0.5  and  r\2  =  1.0  in  0(x)  where  the 
response  turns  unstable  for  the  non-observable  case.  In  subsequent  simulations, 
rj2  is  set  to  be  2.0.  Figure  11  summarizes  the  response  of  Y  and  Yp  for  the 
observable  and  non-observable,  with  n.2  =  2.  The  top  two  graphs  are  for  the 
observable  case,  and  the  bottom  two  describe  the  non-observable  case. 

5.  Figure  13:  Variance  of  the  Sliding  Plane  Switching 

Since  sliding  mode  is  a  high-speed  switching,  or  toggling  method  or 
control,  the  speed  at  which  the  switching  is  to  be  accomplished  must  be  defined. 
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Figure  7.  Observable  Tanker,  T|2=0.5 
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Figure  8.  Observable  Tankers,  T|2=2.0. 
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Figure  9.  Non-Observable  Tanker,  iy=0.5 
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Figure  10.  Non-Observable  Tanker,  rj2=2.0. 
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Figure  11.  Observable,  Non-Observable  Comparison 
of  Y,  Y_  for  n2=2.0. 
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This  switching  speed  can  be  defined  in  terms  of  an  angle.  Note  in  Figure  12  that 
as  the  angle  $  approaches  90  degrees,  the  system  is  being  asked  to  change 
instantaneously.  Another  way  of  saying  this  is  that  the  system  is  ordering  the 
towline  attachment  point  to  move  instantaneously  from  port  to  starboard.  <}>=90 
degrees  is  impossible;  the  designer  has  to  allow  the  system  time  to  react. 

Figure  13  shows  the  observable  and  non-observable  tanker  response  of 
Y,  Yp  for  0=30,  <J)=45,  and  0=60.  For  the  observable  case,  overshoot,  the  amount  of 
oscillation,  and  settling  time  increase  as  (J)  increases.  For  the  non-observable 
tanker,  the  response  destabilizes  between  <J)=45  and  (J)=60.  Clearly  there  is  a  limit 
to  the  switching  speed  of  the  system.  (J)=45  was  used  for  the  standard  tanker. 

6.      Figure  14:  Variance  of  LQR  Weighting 

Figure  14  demonstrates  the  effect  of  varying  the  control  effort  in  the 
linear  quadratic  regulator.  As  a  is  increased  from  0.005  to  0.035,  more  effort  is 
used  to  control  the  towed  vessel  motion.  In  Figure  14,  Y  and  Yp  responses  are 
plotted  for  four  different  values  of  a.  In  all  four  graphs,  the  tanker  is  observable. 
Note  that  as  a  increases,  the  overshoot  in  Yp  decreases,  and  the  settling  time 
remains  about  the  same  as  for  the  standard  tanker.  This  makes  sense;  more  effort 
is  being  expended  to  control  Yp  through  Y^.  The  response  of  Y  remains  essentially 
unchanged. 
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Figure  12.  The  Sliding  Plane. 
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Figure  13.  Y,  Y   Variance  When  Changing  Switching  Angle  8. 
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Figure  14.  Tanker  Response  for  Y,  Yp  as  a  Varies. 
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7.  Figure  15:  Observable  Tanker,  Variance  of  Tp 

Figure  15  demonstrates  the  effect  of  varying  the  time  allotted  for  Y  to 
match  Y^  for  the  standard  tanker.  The  T_  used  for  the  standard  case  was  0.5. 

pc  p 

Varying  Tp  does  not  significantly  affect  the  behavior  unless  the  change  is  in  order 
of  magnitude.  This  result  has  favorable  implications;  if  changing  Tp  doesn't 
radically  affect  the  towed  vessel  response,  a  smaller  motor  can  be  used  to  drive 
the  device  moving  the  towline  attachment  point,  Yp. 

8.  Figure  16:  Variance  of  Towline  Length 

Figure  16  illustrates  how  changing  the  towline  length  will  change  the 
response  of  the  tanker,  both  in  the  observable  and  non-observable  case.  For  this 
simulation,  the  non-dimensional  towline  length,  LL,  was  lengthened  and 
shortened  by  twenty  percent.  If  the  towline  is  shortened,  the  overshoot  and 
oscillation  in  Y  and  Yp  will  increase  for  both  observable  and  non-observable  cases; 
the  reverse  is  true  is  the  towline  is  lengthened.  Settling  time  will  remain  the  same. 

9.  Figure  17:  Variance  of  Towline  Tension 

Figure  17  shows  how  changing  the  non-dimensional  towline  tension 
will  change  tanker  response.  If  the  tension  in  decreased,  the  overshoot  remains 
about  the  same;  oscillation  and  settling  time  in  Y  and  Yp  increase  for  both  the 
observable  and  non-observable  cases.  The  reverse  is  true  as  the  tension  is 
increased.  A  20  percent  change  in  towline  does  not  cause  a  radical  change  in 
system  response. 
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Figure  17.  Observable,  Non-Observable  Tanker:  Variance  of  T. 
T=.0008  T=.0010 T=.0012 
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10.    Non-Linear  (Large  Angle)  Motion 

Recall  from  page  6  that  the  horizontal  towline  force  associated  in  the 
sway  direction  is 

T  sin  ( y  +  \\f ) 
and  the  moment  due  to  the  towline  causing  the  vessel  to  yaw  is: 

-Tx  (  sin  ( y  +  y ) )  -Ty  cos  ( y  +  \\f ) 

Previously  the  assumption  was  made  that  the  range  of  motion  for  y  and  \|/  was 
less  than  30  degrees.  This  simulation  lifts  this  restriction;  the  sine  and  cosine 
terms  from  the  towline  tension  term  are  left  in  when  running  the  simulation. 
Figure  18  plots  the  standard  tanker  response  for  both  linear  and  non-linear 
tensions.  Note  that  the  tanker  response  remains  stable;  the  sliding  mode  control 
law  is  robust  enough  to  handle  non-linearities  in  the  towline  tensions.  The  non- 
linear Fortran  simulation  program  is  included  in  the  Appendix. 
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Figure  18.  Non-Linear,  Linear  Tankers  With  Standard  Parameters. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.      CONCLUSIONS 

This  thesis  explored  the  possibility  of  athwartships  movement  of  the  towline 
attachment  point  as  a  means  of  improving  the  stability  of  towed  vessels.  Newell 
demonstrated  that  use  of  a  full  order  observer  will  work  only  when  the  towed 
vessel  trajectory  is  close  to  that  of  the  towing  vessel  [Ref.  5].  Sliding  mode  control, 
in  conjunction  with  the  Linear  Quadratic  Regulator  cost  function,  provides  a 
much  more  robust  method  of  controlling  the  towed  vessel  motion.  The  sliding 
mode  control  law  can  ensure  stability  of  the  system,  even  if  not  all  the  state 
variables  can  be  observed.  The  only  restrictions  for  the  numerical  applications 
considered  in  this  work  are  that  a  toggling  speed  such  that  (J)  is  less  than  an  angle 
between  45  and  60  degrees  is  used,  and  that  the  toggling  gain  rj2  is  greater  than 
a  number  between  0.5  and  1.0.  The  LQR  function  allows  the  designer  to  choose 
where  his  emphasis  is  placed.  By  choosing  whether  to  place  the  emphasis  on  the 
control  effort,  the  designer  can  minimize  the  size  of  the  motor  required  to  drive 
the  mechanical  apparatus  used  to  move  the  towline  attachment  point.  If  the 
designer  wants  to  minimize  steady  state  error,  he  can  change  the  LQR  function 
to  reflect  that  emphasis.  The  LQR  function  can  also  be  used  to  minimize  specific 
state  variables,  such  as  in  the  case  studied,  where  \j/  was  chosen  to  be  minimized. 
It  should  be  emphasized  that  a  "worst-case"  scenario  was  used  to  demonstrate  the 
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robustness  of  the  system;  a  large,  unwieldy  tanker  was  used,  and  two  of  the  four 
state  variables  were  declared  "not  observable."  Even  in  the  worst  case  the  sliding 
mode  control  law  performed  adequately.  Towline  tension  and  length  were 
increased  and  decreased  by  20  percent,  and  the  system  remained  stable.  The  small 
angle  motion  assumption  was  dropped,  and  yet  the  system  continued  to  perform 
satisfactorily. 

B.       RECOMMENDATIONS 

The  sliding  mode  control  law  should  be  robust  enough  to  cover  most  of  the 
unmodeled  dynamics  and  non-linearities,  but  two  of  the  assumptions  made  at  the 
beginning  of  the  thesis  must  be  investigated  further;  the  massless,  inextensible 
towline,  and  the  no  wind,  wave,  or  current  (disturbances)  assumptions.  In  reality 
the  towline  tension  does  not  remain  constant;  consequently,  the  surge  velocity  u 
of  the  towed  vessel  will  also  vary.  A  method  to  take  these  aspects  of  motion  into 
account  needs  to  be  looked  at;  integral  control  is  one  possibility.  Another  is  to 
determine  a  wave  and  towline  function,  and  add  the  two  functions  to  the  right 
side  of  the  sway  and  yaw  equations  of  motion. 

This  thesis  has  shown  that  it  is  possible  to  stabilize  and  control  the  motion 
of  a  towed  vessel  with  a  sliding  mode  control  law  that  uses  linear  quadratic 
regulator  gains.  Further  analysis  needs  to  be  done  on  the  effects  of  disturbances 
and  a  time-variant  towline  tension  on  the  system. 
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APPENDIX 


TABLE  I.  TOWED  VESSEL  DATA 

Vessel 

Property 

Barge 

Tanker 

Mariner 

LBP,  ft 

191.56 

1066.3 

528 

4 

-0.00136 

-0.0009 

-0.0004444 

m' 

0.170 

0.0181 

0.00888 

y' 

1  V 

-0.01383 

-0.0171 

-0.00912 

y' 

-0.0153 

-0.0261 

-0.01434 

y' 

0 

0 

0 

Y> 

0.00238 

0.00365 

0.00456 

Ni 

0 

0 

0 

N'v 

-0.007285 

-0.0105 

-0.0046 

ii-Ni 

0.00188 

0.00222 

0.00115 

Nf 

-0.00128 

-0.0048 

-0.00296 
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TABLE  2.  NONDIMENSIONAL  TERMS 


=  m  I  mr 

-  u  /  ur 

=  v  /  uc 

=  u  /  (u2c-L) 

=  v/   (ul-L) 

=  t /  (  ut  - L  )       (  time  ) 

=  IJ  (mrL2  ) 

=  \|/  =>]//  (  u0/L) 

=  V|>/  =  V|//   {ul  /  L2  ) 


*g  =  *g  /  1 


yg  =  yg  /  i 

x'  =  XU/  (  mruc  /  L'  ) 

x'  =XJ  mr 


Y'V  =  YV/  (mruc/  L) 

yJ«Sf«Y?/  (mrL) 

N'V=NV/  (mruc) 

Ni=Nj  (  mrL  ) 

tfj  =i<  =A^/  (  mr  ucL) 

N't  =Ny=Ny/   (mrL2  ) 

T'  =  T  /  (  mrul  /  L)      (  Tension  ) 

1'  =  2/  L 

x'F=xp/  L 

y'P=yJ  l 
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MATRIX,  LQR  PROGRAM 

Xudot=- 0.0009; 

M=0.0181; 

Yvdot=-0.0171; 

Yv=-0.0261; 

Yrdot=0.0; 

Yr=0.00365; 

Ydel  =0.0278; 

Nvdot=0.0; 

Nv=-0.0105; 

Nrdotlz=-.00222 

Nr=-0.0048 

Ndel=-0.0139 

inquire  LL 

inquire  Xp 

inquire  T 

U=1.0; 

L=1.0; 

D=NrdotIz*(Yvdot-M)-Nvdot*Yrdot; 

al3=1.0; 

a24=1.0; 

a31=(l/D)*(T/LL)*(NrdotIz-Xp*Yrdot); 

a32=(l/D)*((Yv*U+T*(l+Xp/LL))*NrdotIz-Yrdot*(Nv*U+'rXp*(l+Xp/LL))); 

a33=(l  /D)*(Nv*Yrdot-Yv*NrdotIz); 

a34=(l/D)*((Yvdot*U-Yr)*NrdotIz-Yrdot*(Nvdot'fU-Nr)); 

a41=(l/D)*(T/LU*(Xp*(Yvdot-M)-Nvdot); 

a42=(l/D)*((Yvdot-M)*(Nv*U+T*Xp*(l+Xp/LL))-Nvdot*(Yv*U+T*(l+Xp/LL))); 

a43=(l/D)*(Nvdot*Yv-Nv*(Yvdot-M)); 

a44=(l/D)*(Nvdot*Yr-Nr*Yvdot+(Nr-Nvdot*U)*M); 

b3=(l/D)*T(NrdotIz/LL-Yrdot*(l+Xp/LD); 

b4=(l/Dn>((Yvdot-M)*(l+Xp/LL)-Nvdot/LL); 

A=[0,0^13,0A0/0^24^31/a32/a33,a34;a41,a42/a43,a44]; 

B=[0,O;b3;b4l; 

C=[0,0,1,0]; 

D=[0]; 

Q=[0  0  0  0;0  131.3316  0  0;0  0  0  0;0  0  0  0]; 

inquire  alpha 

R=[(l/alpha)**2]; 

<EIG,K>=REGULATOR(A,B,Q,R); 

INQUIRE  TP 

Al=[-l/TP/0,0/0,0;0,0,0,l,0,-0,0,0,0,l;B3/A31/A32,A33/A34;B4/A41,A42,A43,A44]; 

Bl=[l/TP;0/-0/O;0]; 

S=[1,K]; 

G1=-INV(S*B1)*S»A1; 

G2=-DVJV(S*B1); 

ETA2=G2 

GAtf«J=Gl 

SIGMA=S 
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LINEAR  SIMULATION  PROGRAM 


C 
C 


C 


REAL  KYP,K1  ,K2,K3,K4,L,LL,Nvdot,M,Nv,Nr,Nrdot 

Xudot=-0.0009 

M=0.0181 

Yvdot=-0.0171 

Yv=-0.0261 

Yrdot=0.0 

Yr=0.00365 

Ydel=0.0278 

Nvdot=0.0 

Nv=-0.0105 

Nrdot=-0.00222 

Nr=-0.0048 

Ndel=-0.0139 

U=1.0 

L=1.0 

OPEN  (10,FILE='LFNEAR.DAT',STATUS='OLD') 
OPEN  (1 1,FILE='Y.DAT',STATUS='NEW') 
OPEN(12,FILE='PSI.DAT',STATUS='NEW) 
OPEN(13,FILE='YP.DAT',STATUS='NEW') 
OPEN(14,FILE='YPC.DAT',STATUS='NEW') 

READ  (10,*)  TSIM,DELTAT,IPRNT 
READ  (10,*)  KYP,K1,K2,K3,K4 
READ  (10,*)  SYP,S1,S2,S3,S4 
READ  (10,*)  GBAR 
READ  (10,*)  TP,LL,XP,T 
READ  (10,*)  ETA2,PHI 
READ  (10,*)  Y,PSI,SAT 

ISIM=TSIM/DELTAT 

PI=4.0*ATAN(1.0) 

PHI=PH1*PI/ 180.0 

PSI=PSI*PI/ 180.0 

SATP=  SAT 

SATM=-SAT 

XI  =Y 

X2=PSI 

X3=0.0 

X4=0.0 

YP=0.0 

YPC=0.0 

D=Nrdot*(Yvdot-M)-Nvdot*Yrdot 
a!3=1.0 
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a24=1.0 

a31=(l/D)*(T/LL)*(Nrdot-Xp*Yrdot) 

a32=(l/D)*((Yv*U+T*(l+Xp/LL))*Nrdot-Yrdot*(Nv*U+TXp*(l+Xp/LL))) 
a33=(l  /D)*(Nv*Yrdot-Yv*Nrdot) 
a34=(l/D)*((Yvdot*U-Yr)'Nrdot-Yrdot*(Nvdot*U-Nr)) 
a41=(l/D)*(T/LL)*(Xp*(Yvdot-M)-Nvdot) 
a42=(l/D)*((Yvdot-M)*(Nv*U+T*Xp*(l+Xp/LU) 
&     -Nvdot*(Yv*U+T*(l+Xp/LL))) 
a43=(l  /D)*(Nvdot*Yv-Nv*(Yvdot-M)) 
a44=(l/D)»(Nvdot*Yr-Nr*Yvdot+(Nr-Nvdot*U)*M) 
b3=(l/D)*T*(Nrdot/LL-Yrdot*(l+Xp/LL)) 
b4=(l/D)*T((Yvdot-M)*(l+Xp/LL)-Nvdot/LL) 
C 

C       SIMULATION  BEGINS 
C 

DO  1  I=1,ISIM 
YPDOT  =-YP/TP+YPC/TP 
X1DOT  -  X3 
X2DOT  =  X4 

X3LXDT  =  B3*YP  +  A31*X1  +  A32*X2  +  A33*X3  +  A34»X4 
X4LX)T  =  B4*YP  +  A41*X1  +  A42*X2  +  A43*X3  +  A44*X4 
C 

YP  =  YP  +  YPDOT*DELTAT 
XI  =  XI  +  X1DOTDELTAT 
X2  =  X2  +  X2DOT*DELTAT 
X3  =  X3  +  X3DOTDELTAT 
X4  =  X4  +  X4DOTDELTAT 
C 

C  CONTROL  LAW 

C 

SIGMA  =  SYP*YP  +  S1*X1  +  S2*X2  +  S3*X3  +S4*X4 
SATSGN=  SIGMA/TAN(PHI) 
IF  (SATSGN.GE.(+1.0))  SATSGN=+1.0 
IF  (SATSGN.LE.(-l.O))  SATSGN=-1.0 
UHAT   -  KYP*YP  +  K1»X1  +  K2*X2  +  K3*X3  +  K4*X4 
UBAR   =  ETA2*GBAR*SATSGN 
YPC     =  UHAT+UBAR 
IF  (YPC.GE.SATP)  YPC=SATP 
IF  (YPC.LE.SATM)  YPC=SATM 
C 

Y   =X1 
PSI-X2 
C 

C  PRINT  RESULTS 

C 


H+l 

IF  (J.NE.IPRNT)  GO  TO  1 
TIME=I*DELTAT 
WRITE  (V)  TIME,Y 
WRITE  (11,*)  TIME,Y 
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WRITE  (12,*)  TIME,PSI*180.0/PI 

WRITE  (13,*)  TIME,YP 

WRITE  (14,*)  TIME,YPC 

J=0 

CONTINUE 
STOP 
END 
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NON-LINEAR  SIMULATION  PROGRAM 


C       NONLINEAR  SIMULATION  -  TANKER 

C       SURGE  NOT  INCLUDED 

C 

RE AL  KYP,K1  ,K2/K3,K4,L,LL,NVDOTNV,IZ,NRDOT,NR,NVRR 
C 

XUDOT=-0.0009 

M      =  0.0181 

YVDOT=-0.0171 

YV     =-0.0261 

YRDOT=  0.0 

YR     =  0.00365 

NVDOT=  0.0 

NV     =-0.0105 

IZ     =0.0 

NRDOT=-0.0022 

NR     =-0.0048 

YVRR  =-0.045 

NVRR  =  0.0061 

UTOW  =  10.0 

L      =  1066.3 

RHO  =  1.9905 

SB     =  680625.0 

P      =  9.78 

Q      =  1.93 

SB    =  SB/(0.5*RHO*L*L*UTOW*UTOW) 
C 

OPEN(10,F1LE='LINEAR.DAT',STATUS='OLD') 

OPEN  (1 1,FILE='Y.RES',STATUS='NEW) 

OPEN  (12,FILE='PSI.RES',STATUS= 'NEW') 

OPEN  (13,nLE='YP.RES',STATUS='NEW') 

OPEN(14,FILE='YPC.RES',STATUS='NEW') 

OPEN  (15,FILE='NONLINEARl. RES' ,STATUS='NEW) 
C 

READ  (10,*)  TSIM,DELTAT,IPRNT 

READ  (10,*)  KYP,K1,K2,K3,K4 

READ  (10,*)  SYP,S1,S2,S3,S4 

READ  (10,*)  GBAR 

READ  (10,*)  TP,LL,XP,T 

READ  (10,*)  ETA2,PHI 

READ  (10,*)  Y,PSI,SAT 

ISIM  =  TSIM/DELTAT 
PI     =  4.0*ATAN(1.0) 
PHI   =  PHI*PI/180.0 
PSI   =  PSI*PI/ 180.0 
SATP  =  SAT 
SATM  =-SAT 


C 
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X  =  LL-XP 
U  =  1.0 
V  =  0.0 
R  =0.0 
YP  =  0.0 
YPC   =0.0 

c 

DEN=(M-WDOT)*(LZ-NRDOT)-YRDOT*NVDOT 
C 

C       SIMULATION  BEGINS 
C 

DO  1  I=1,ISIM 
C 

F2=YV»V+0.5*YVRR*V*R*R+(YR-M)*R 

F3=NV*V+0.5*NVRR*V*R*R+NR*R 
C 

SING=(Y+XP*SIN(PSI)+YP*COS(PSI))/LL 

GAMMA=ASIN(SING) 
C 

YPDOT  =-YP/TP+YPC/TP 

VDOT   =(F2-T*SIN(GAMMA+PSI))/(M-YVDOT) 

RDOT=(F3-T*XP*SIN(GAMMA+PSI)-T*YP*COS(GAMMA+PSI))/(IZ-NRDOT) 

YDOT  =U*SIN(PSI)+V*COS(PSI) 

PSIDOT=R 
C 

YP   =  YP   +  YPDOT  *DELTAT 

U     =  U     +  UDOT  *DELTAT 

V  =  V  +  VDOT  *DELTAT 
R  =  R  +  RDOT  *DELTAT 
X     =  X     +  XDOT  *DELTAT 

Y  =  Y  +  YDOT  *DELTAT 
PSI  =  PSI  +  PSIDOT*DELTAT 

C 

C  CONTROL  LAW 

C 

SIGMA  =  SYP*YP  +  S1*Y  +  S2*PSI  +  S3*YDOT  +  S4*PSIDOT 

SATSGN=  SIGMA/TAN(PHI) 

IF  (SATSGN.GE.(+1.0))  SATSGN=+1.0 

IF  (SATSGN.LE.(-l.O))  SATSGN=-1.0 

UHAT  =  KYP*YP  +  K1*Y  +  K2*PSI  +  K3*YDOT  +  K4*PSIDOT 

UBAR  =  ETA2*GBAR*SATSGN 

YPC     =  UHAT+UBAR 

IF  (YPC.GE.SATP)  YPC=SATP 

IF  (YPC.LE.SATM)  YPC=SATM 
C 

C         PRINT  RESULTS 
C 

J=J+1 

IF  (J.NE.IPRNT)  GO  TO  1 

TIME=I*DELTAT 
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WRITE  (*,*)  TIME,Y 

WRITE  (11,*)  TIME,Y 

WRITE  (12,*)  TIME,PSI*180.0/PI 

WRITE  (13,*)  TIME,YP 

WRITE  (14,*)  TIME,YPC 

WRITE  (15,*)  TIME,Y,PSI*180.0/PI,YP,YPC 

J=0 
CONTPNUE 
STOP 
END 


54 


LIST  OF  REFERENCES 

1.  Bernitsas,  M.M.,  and  Kekridis,  N.S.,  "Simulation  and  Stability  of  Ship 
Towing,"  International  Shipbuilding  Progress,  v.  32,  1985. 

2.  Papoulias,  F.A.,  "A  Qualitative  and  Quantitative  Study  of  Steady-State 
Response  of  Towed  Floating  Bodies,"  Dynamics  and  Stability  of  Systems,  v.  3, 
1988. 

3.  Kolthoff,  D.L.,  "Parametric  Study  of  the  Dynamic  Stability  of  Towed  Ships," 
Master's  Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  December,  1989. 

4.  Krikelis,  N.J.,  and  Kavouras,  D.,  "Dynamic  Performance  of  Towed  Vessels 
Employing  Nonlinear  Rudder  Control  Under  Adverse  Weather  Conditions," 
International  Shipbuilding  Progress,  v.  33,  1986. 

5.  Newell,  J.B.,  "Automatic  Control  of  Straightline  Motions  of  Towed  Vessels," 
Master's  Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  March  1990. 

6.  Friedland,  B.,  Control  System  Design,  McGraw-Hill,  Inc.,  1986. 


55 


INITIAL  DISTRIBUTION  LIST 


No.  of  Copies 

1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library,  Code  52  2 
Naval  Postgraduate  School 

Monterey,  CA   93943-5002 

3.  LT  A.K.  Samora  1 
255  Bethany  Dr. 

Scotts  Valley,  CA  95066 

4.  Naval  Engineering  Curriculum  Office,  Code  ME  1 
Naval  Postgraduate  School 

Monterey,  CA   93943 

5.  F.A.  Papoulias,  Code  ME /Pa  3 
Naval  Postgraduate  School 

Monterey,  CA   93943 

6.  Chairman  1 
Department  of  Mechanical  Engineering,  Code  ME 

Naval  Postgraduate  School 
Monterey,  CA  93943 


56 


Thesis 

S1555  Samora 

c.l      Sliding  mode  control  of 
motions  of  towed  ships. 


^ 


